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ABSTRACT 


A method is reviewed for computing the aerodynamic influence coefficients 
{AlCs) for slender bodies. The method is based on the unsteady slender- 
body theory by Miles and its extension to obtain the AICs by Rodden and 
Revell. 


The simplicity of the slender-body theory permits the definition of a number 
of sets of AlCs for use in transient analysis. The influence coefficients 
relating the transient aerodynamic forces to the body deflections and their 


first two derivatives are defined by the following relation: 
a Saeey "2 ,..2 
{F(t)} = (gs/arc, 14ht + [c, line/v} + [c, J {ne lV ) 


The matrices [c, .] ; [c, 4}. and [C, 5] are seen to be steady, damping, and 
inertial AlCs, respectively. The oscillatory AIlCs are defined by 


2 


{r} = pw bés[ C,] {ht 


and are related to the above definitions through 
ak“ (és/s)C,] =[C,_] + ik (é/b \LC, 4) - K2(e/b_)*LC,.] 
r h hs r r’- hd r r hi 


The Aerospace IBM 7090 Computer Program No. HM15 provides the AICs 
in printed and optional punched-card output formats. The program capacity 
is 50 control points and, in the oscillatory case, 50 values of reduced 


velocity. 
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SECTION I 
FORMULATION OF PROBLEM 


A; Introduction 


The simplest theory available for the estimation of unsteady aerodynamic 
loads on slender bodies is the momentum theory of Munk! as extended to the 


Die The limitations of the so-called slender-body 


unsteady case by Miles. 
theory have also been summarized by Miles (Ref. 4, Table 2, p. 161). These 
are: (1) the fineness ratio must be much less than unity; (2) the Mach number 
must be much less than the reciprocal of the fineness ratio; and (3) the motion 
must be slowly varying, or, more specifically, in the oscillatory case the re- 
duced frequency must be of order unity. Hence, we see that the slender-body 
theory can be useful in the supersonic flight regime provided that the body is 


sufficiently slender and the motion is not violent. 


The present formulation is based on the derivation of the oscillatory 
aerodynamic influence coefficients (AICs) from slender-body theory given by 
5 i 
Rodden and Revell. 


tension of Ref. 5 to.obtain a number of sets of AICs for use in transient 


The simplicity of slender-body theory permits the ex- 


analysis, and this extension is made in the present treatment. 


B. Sign Convention 


The standard NASA stability axis system is used throughout. The positive 
directions of x, y, and z are forward, starboard, and downward, respectively. 
Positive rotations are given by applying the right-hand rule to the coordinate 
directions. The vehicle is assumed to be moving in the positive x-direction; 
i.e., the relative wind is in the negative x-direction. The force sign con- 


vention is the same as the coordinate convention. 


Cc. Derivation of Equations 
The derivation of equations for the oscillatory case is given in Ref. 5 
(pp. 60-72). However, Ref. 5 containsa mistake in sign that leads to the 


incorrect sign on the damping terms, so it becomes necessary to repeat the 


derivation here with the correction. We shall present the derivation for the 
oscillatory case since, by properly identifying the various terms, we may 
indicate the results for the transient case. We define the transient influence 
coefficients to relate the aerodynamic forces to the body deflections and their 


first two derivatives by the following 
{F(t)} = (qs/2)({c, 1h} + [c, {he/v} + Ic, Wie v"4) ; (1) 


The matrices [Cel [C, al and [C, are seen to be steady, damping, and 
inertial AICs, respectively. The oscillatory AICs are defined by 


{F} = po“b“s[C, ]{h} . (2) 


Equations (1) and (2) must be identical in the oscillatory case, so the following 


relationship must exist among the various AICs 
ak2(és/s)(C,] = [C, _] + ik (é/b_)[C, ,} - k2(é/b_)* Ic, |] (3) 
r h hs r r'- hd r r hi : 


The equations given by Bisplinghoff, Ashley, and Hialiniane provide a 
convenient basis for deriving the AIlCs of a slender body. The vertical force 
acting per unit length of the body is the reaction to the substantial rate of 


change of the momentum of the virtual mass per unit length of the body 


#--2(B) : 


ae Verte (5) 


since the relative wind is moving in the negative x-direction (see Fig. 1). 


Vv 


RELATIVE 
WIND 


Fig. 1. Slender-body geometry for AICs. 


The momentum of the virtual mass per unit length is found from the local 


cross-sectional area and the downwash 


& = pSlx)w(x, t) (6) 


where the downwash is the substantial derivation of the deflection 


w(x, t) = De (7a) 
Z oh , Oh 
=-VERt (7b) 
If we assume harmonic motion and substitute Eqs. (5), (6), and (7b) into 
Eq. (4), we obtain 
dF yd dh, . af dh, . 
ae pV $ |S (-V dx + ih) | - iwpS(x)(-V as + iwh) . (8) 


To obtain the force on a specified length of the body, it is necessary to inte- 


grate Eq. (8) over that length. We consider the body to be divided into a 


“number of sections, not necessarily of equal length, as shown in Fig. 1. For 


the n'th section, the control point is taken at the midpoint of its length An: ¥ 


The aft end of the section is located at x 1/2 ae ieee a,/2 and has the cross- 
section area Se 1/2' the forward end of the section is located at 


ree + A,/2 and has the cross-section area S Carrying out the 


x : 
n+1/2 n+1/2 
integration of Eq. (8) for the n'th section yields the following 
x 
FO - dx dx (9a) 
*n-1/2 


nt 


pv [ s(x) (-vs2 ion) D422 : wf eee (vee + iwh ) dx (9b) 


=e [-V'Spoiy2 hh41/2 7 Sn-1/2 ®n-n/2)t (8 a41 72 Pag /2 ” $,-1/24n-1/2)| 
x 
: wn [ mt 1/2 a y(-Vh! + wh) dx (9c) 
*n-1/2 


We resort to Lagrangian interpolation for the evaluation of the terms in 


Eq. (9c). For the first two terms, we choose parabolic interpolation 


n n n n+] 
(x -x  ,)(x- x) 
ers | x es x - x (10) 
+1 n-l'*' ntl n 


a 

“The input to the program must be the coordinates of the endpoints (aft 
and forward) of each section, from which the control point locations are cal- 
culated by averaging. With a given total body length, it is not possible to 
specify the control point locations a priori. 


Be 


For evaluation of the integral, we use only linear interpolation 
h=h + hi(x - x) (11) 
and 
Ste) = ANS, 1/2 pary2 ~) * Spa ipal® pera} 2) 


The three-point interpolation leads to three elements in each row of the 


AIC matrix. Therefore, from the definition of Eq. (2), 


2,2 
Fa eilica b(n nel nel - Cin hn bi Cry, ntl Ont)! ; (13) 


By identifying Eq. (13) and Eq. (9c), after evaluating the integral of Eq. (9c) 
by means of Eqs. (11) and (12), we obtain the following relation for the AICs 


= 2 
Ohare wed ss Cie nd * Con net nel ~ (1/25) (Sasi bhi9a> Sena tae 


+ i(k /D S472 2aea/2 7 Sa-1/22n-1/2 + AV a) 


(E/E) [(On/12) Saay2 - Speiy2 a BR OV,}} 6 04) 


where AV. is the volume of the n'th section. 


If we evaluate the deflections in Eq. (10) at stations x and x 


n-1/2 
and then, by differentiating Eq. (10), evaluate the slopes at stations x 


n+1/2’ 


n-1/2’ 
then we may place the right-hand side of Eq. (14) in terms of the 


h_, andh 
n n 


Xn’ nt) /2? 


control point deflections hooey +L If we carry out this substitution, 


the AICs are found by identifying the coefficients of the control point deflec- 
tions on both sides of Eq. (14). The following results are obtained 


= (unten) 3 
esl _ (1/kesD. 2541/2! 4nt1 7 4.) 25,,-1/2'34y - Ant) - i(k /b,) 


x |- Sn41/24n Anti ~ Sy-1/2 4nl@4n + Sngy) - 2AV AIA, + O07) 
- (RE/PE) (45/6) Saary2 - Saeas2nt Ana)| ite) 
where 
Diy = (4y_) + AMAL) + 24) + AL) ; (16) 
ign : (resp, ) {2841 ;2l4q 1 bi 2A, - Ant) . 2S-1/2!4n-1 7 2A, ~ Any)) 
+ i(k, /b,)|s nt]/2 Ans pOn-1 + 24,) - Sins 1/24 _j(24, +4 11) 


. 2AV (Ant az 4,_,)| 7 (x21?) [(a2/6) (S42 ~Sigr/2) Ay hl An-1) -p, av, 


(17) 


where 


+A VA + CD) ; (18) 


2 2 ){ ; : aha 
ee te | (1/kfsD. 2S-1/2'4n-1 7 4.) ~ 28 41/2'9nert 3A4,) Z i(k,/b,) 


a [Ssasie ey ss 2a.) m SE Seo neo + 2AV (Ay + 4,) 


2 75.2 2 
- (Kr/b)) (05/6) ety? Page vel OR! , (19) 


where 


ae (A-4 + 2A, + Anal, tA (20) 


n+ +) 


The above expressions are applicable for all intermediate sections of the body, 
sections which are centrally located as far as the interpolation and differen- 


tiation are concerned. The exceptions are the first and N'th sections. The 


counterpart of Eq. (14) for the first section is 
= 2 1 _ t . 
Cnr. t Cnr, 2%2 * Sn, 323 = (1/2) [si j2 4 yo S32 hgy2 + Uk, /b,) 


X (S83 72h372-Syp2hysothyAv)) - (xi /bF )I(aq/ 12) (S72 $372) hy -h, av, || 


(21) 


Carrying out the evaluation of the appropriate deflections and slopes from 
Eq. (10) in terms of the first three control point deflections leads to the fol- 


lowing coefficients for the first row of the AIC matrix 


- (unden,)| : 
Chi, 1° (i/K2sp, 283 72(34, + A3) - 2S) 72(44, +3A,+ A3) + i(k ./b,) 
x [55,2 49(2a,+4;) ~Syp2(2A, + ANZA, + 2A, +43) - 2AV (2A, + 3A, +4,)| 


(x22) [ (4/6) (S372 - Syp2(2@a, + 342 + 43) - p,av,|| % Wee) 


where 


EDs 


tthe, ag: , 
Orn he (1/K2sp, 284 )2(A, ~ 2A, A) + 28 ,2(3A, +20, +A,) + ilk, /b,) 


x [s3).4,(24, + A,)+ 51724 ,(24, + 2A, + A) + 2AV (A, + 2A, + 43) 


- (27) (44/6) (S,)2 - S3,.MA, +20, + al, (23) 


where 
D,=D_ withn =2 : 
2 n 
: 2 )| : 
Cats (1/K2sp, 283 72(49 - A) - 28, ,9(3A, +A.) + i(k, />,)| 8455 Ay As 


: $1/2.4(24, +4,) -2av,(a, +4,)] : (x5 /F) (A576) (S39 -8, pM, + AQ) 


(24) 
where 
D, = Daisy with n = 2 
Similarly, the counterpart of Eq. (14) for the N'th section is 
Cc hoes h ,+C ie.2 (1/K2s) {s hi 
hN,N-2 N-2 hN,N-1 N-1] hN,N'N r N-1/2 “N-1/2 
a ' ; = t 
Su4i/2 N+1/2 + Hc MSy 572 "N12 Sy-1/2%N-1/2 t hy Vy) 
. (K2/02 (ax,/12) (S as yhi - hav |} (25) 
ror N N-1/2 N+1/2’ N NON ' : 


=H 


and we obtain the following coefficients for the last(N'th) row of the AIC 


matrix 
C = (1/K2sp,,_ 4) {2s (Acc ahs 2s (BA + Aggy) 
hN, N-2 oo oN 2) [N12 Ned ON N+1/2 N-1 
+ ilk, LS 4 1/2 ADAG A Ae aS erry AA 1M eAV A Pale 4} 
2,,2\/.2 
: (5/2) (an /6) (Sye172 7 Sngis2i4n t Andy? (26) 
where 
Dy-2 = Py-y Withn=N- : 
(A a & 
ChN,N-1 =(1/K2 sDy 4) (28a ps OAn $2An yy tAy_2) + Ay y2ldy > 24n 21 > Ang) 
Wetoichagel Ay (2Ay + 2An yt An a) - Syiys24n(@4nig t Aya) 
- 2A), (Ay, + 2A. | - (x2/02) (2/6) (S -S YA. + 2A +A ) 
4n- 2 r r/\°N N+1/2) “N-1/2"°°N N-1 N-2 
(27) 
where 
Dy.) 7D, withn=N- 1 : 


2 
Cun, N = (1/2 sDyp) [289 _ 1 72(3 Aq 1 + Anya) ~ 28yq4 1/244 + 3AN 1 + Ana) 


+ i(Ic, /b,) [Swany2l2On + Any. An + 24Q 1 + Ag) - Sy y2 An (24-1 + An. 2) 


+ 2AV (ZO + 34y yt Ay.) ‘ («E/E (an /6) (Sy a/2 ~ Syi12) 


(2A, + 3An y+ An2) - Dy AV | 


(28) 


where 


with n=N- 1 


io) 
u 
o) 


N “ntl 


To illustrate the assembly of the coefficients into the AIC matrix, we 


show the format below for a slender body having five degrees of freedom. 


Chick ~~ SRE CHS 0 
yee Wa ey Chad Q Q 
[o.] = Ona Cease.  “Sygce e (29) 
9 Ch4,3 Cn4,4 Sn, 
° ? “ns,3 Sns,4 Ons, 5 


If we formally write each element of the oscillatory AIC matrix in the 


form 
: 2 2,2 
Cc a (usZep,)[A,, + i(k, /b,) By - (Ki/b2) ¢,,| (30) 


then the elements of the transient AIC matrices may be identified by compari- 
son with Eq. (3). We see that the elements of the steady AIC matrix are 
given by 


hs. 
qJ 


Cc _ = 2(t/8) A;,/D, : (31) 


the elements of the damping AIC matrix are given by 


Ona, = (2/S) B,,/D; ; (32) 


-10- 


and, finally, the elements of the inertial AIC matrix are given by 


eo = (2/&S) C,./D, : (33) 


Before concluding this discussion, we note that, according to Miles 
(Ref 4, p.169) the slender-body theory presented above may be applied to 
finned vehicles if an effective cross-section area in the finned region, whose 


geometry is shown in Fig. 2, is taken as 


2 


Ss, * also Re + RA s*) (34) 


Fig. 2. Cross section of finned 
region of slender body. 
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SECTION II 
GENERAL DESCRIPTION OF INPUT 


A. Units 


Since all dimensional input is geometrical and the aerodynamic matrix 
is dimensionless, only a consistent set of length units is necessary: inches 


or feet. 
B. Classes of Numerical Data and Limitations 


The data required by the program are control indicators, geometry, 


and a set of reduced velocities. The example problem illustrates their use. 


L. Example Problem 


We consider the five-segment body shown in Fig. 3 computing the 
oscillatory case for the reduced velocities (1/k)) of 2.0 and 6.0, the transient 


case (steady, damping, and inertial AICs), and the steady case (steady AICs only). 


$:10.0 $=12.0 $=14.0 $=1!,0 S= 


Fig. 3. Five-segment body. 


The geometrical data required as input to the program are listed 


below. 
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Segment No. x (ft) S (sq ft) 


1 (aft) 1 10 
1 (forward), 2 (aft) 3 12 
2 (forward), 3 (aft) 6 14 
3 (forward), 4 (aft) 10 ll 
4 (forward), 5 (aft) 12 8 
5 (forward) 15 0 

c = 10.0 ft 

b = 5.0 ft 

r 

s = 20.0 ft 

S = 200.0 sq ft 


2. Program Restrictions and Options 


a. The maximum number of segments into which a body may be 


subdivided must be < 50. 

b. The maximum number of 1/k 's for one data deck must be 
< 50. 

c. If it is desired to compute the steady case, a zero value of 


V/k must be used, to compute the transient case a negative value of V/k 


must be used. 
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SECTION II 
DATA DECK SETUP 


Ay Loading Order 


Input data decks punched from keypunch forms are loaded behind column 


binary deck HM15. Any number of complete decks may be stacked. The data 


for each deck must be in the following order: 


(1) 
(2) 


(6) 


Heading card 

IX, KVBRW, IPUNCH 
CBAR, BR, S, SREF 
X(I) series 

LOCAL AREA(I) series 


(1/k,), series 


B. Input Data Description 


(1) 


(2) 


The heading card may contain any information desired in Columns 


2 through 72. Column 1 is always left blank. 
Control card (FORMAT 1814) 


(a) IX = the number of control points; also the number 


of body segments, < 50 


(b) KVBRW = the number of reduced velocities to be listed 
in data item (6), < 50 


(c) IPUNCH 


0 or blank if the computed matrices are to be 


punched in cards; 


# 0 if no punched output is desired. 
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(3) Single parameters card (FORMAT 6E12. 8) 
(a) CBAR = c, reference chord 
(b) BR = bi. reference semichord 
(c) S = s, reference semispan 
(d) SREF = 5S, reference area 


The constants are tabulated in the order in which they are defined. 
CBAR and SREF are used in the steady and transient cases; BR 
and S (semispan) are necessary for computing oscillatory cases. 
When only one pair of the constants is needed, the two respective 


card fields for the other pair may contain zero or be left blank. 
(4)  X{(I) Segment coordinate series, I <5] (FORMAT 6E12. 8) 


The x coordinates (x and x ) used in this series locate 
n-1/2 n+1/2 
the aft and forward end of each body segment. The number of terms in 
the series is one more than the number of control points (I= 1, IX #1). 
Begin the series with %1/2 (aft end) of the aft body segment and list . 
the consecutive coordinates through Xl /2 (forward end) of the foremost 


body segment (nose). 
(5) LOCAL AREA (I) series (FORMAT 6E12. 8) 


These are the local cross-section areas (5 and S ) at the 
n-1/2 nt+1/2 
aft and forward end of each body segment. The number of terms in this 
series is the same as in X(I). The areas are listed in the same order 
as the coordinates; begin with nad 7 for the aft body segment and list 


the consecutive areas through S for the foremost body segment. 


nt1/2 
(6) (1/k,), series (FORMAT 6E12. 8) 


This series consists of the reference reduced velocities for the 
oscillatory case and the codes for obtaining the steady and transient 


cases. To obtain the steady matrix [Cc 6] input V/k. = 0.0 and for the 


4 
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transient case (Lc, JJ, [c, 4]: and Le, .)) use 1/k = any negative number. 


The number of 1/k,'s in the series is set by KVBRW in the control card. 


NOTE: Each new series starts ona new line (card). 


C. Example Keypunch Forms 


Example keypunch forms are given on the following pages. Columns 73 
through 80 are reserved for data deck identification. This space may be used 
in any fashion; however, it is suggested that the last three columns be used 
for sequencing. Inthe example that follows, only the sequenced cards (lines) 
are to be used inthe sample data deck; the lines with blank Columns 73 through 


80 are for explanation of the input. 
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eto HADAAAAAAAVAANANAANORAAAAATAANAAANOTAONATAANTNT CLL eet TT] 
tN HH if 1 TTT ee 
anit AT Ht Tit i EEE Tse 
EEA } PET 
EET TT Te rt Ee ll Pele CTT bel TTT 
settaiteh TET MTH Ata 
effet sbaettotel LULL | Isteetetote LUT | fel lfetett TET dtetelobd FTE TU Pte CETTE TT Tolle ETT TTT el fl 
UTE ETT tsetse beard eeebebasb id setselcpi] ETT TTT ETT] 
VUOTOUNEEVONOWEWONONONONNQANAOAOSONOTANONOTOVOVOVNTENNOROTOTATOTOTOVOTOGNTEY 
Fe tReet cope 
FLEA ll HARA! tHE PSHE eT 
Ht Ht HH AEE ETE LAH 

ATE Ht He HA I I HATH HE rad 
of fol elt il Neer Ae Phe 
i tt i Ht +} ee PET 44 ble phd pret 

AT EE aucnar al 
PATTI a oe a 
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ti: 


o 
‘eH 
Hida! 


Lee 


a me 


HH tt TERRACE 
Se 


ieee etait 


— . i 


He 


LO 


SECTION IV 
PROGRAM OUTPUT 


A. Printed Output 
l. All input data. 


2. The segment lengths | DELTA X (I)] and the controi point 


coordinates XN (I), (4, and x, in Section I.C). 


3. Each AIC matrix preceded by the associated 1/k., and the matrix 
size. 
4. Sequencing number of the first and last punched cards (output) 


for each AIC matrix. 


5. Example problem printed output is given on the following pages. 
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B. Punched Output 


1. A deck of punched cards (output) from this program is suitable as 


an input deck to other programs requiring AICs. 


2. All punched output is sequenced in order on Columns 73 through 80 


starting with HM150001. The data is punched in the following order: 


a. Card 1 contains (1/k_j: (FORMAT 6E12. 8). 

b. Card 2 contains IX, the size (number of control points) of the 
AIC matrix (FORMAT 1814). 

c. The AIC matrix punched in column binary form and its TRA 


card make up the remainder of the punched output for (1/k)). 


3, The order of Statement 2 can be repeated for all reduced velocities 
per input deck if we note that in the transient case, the three matrices [c, J]. 
[C, 4 ’ and [c..] , 
in the order listed with a TRA card after each matrix. 


which must be considered in Statement 2.c, are punched 


4. Each AIC matrix is punched by columns. Column | starts in 


Orgin 1 and Column 2 in Location (1 + matrix size). 


5. The oscillatory matrices are punched in the order Column 1 (real), 
Column 1] (imaginary); Column 2 (real), Column 2 (imaginary); ...; Column IX 
(real), Column IX (imaginary). Inthe steady, damping, and inertial matrices 


all columns are real and are punched in order. 
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SECTION V 
PROCESSING INFORMATION 


Operation 

STANDARD FORTRAN MONITOR system 
Estimated Machine Time 

T = time in minutes 

IX = number of control points 

KVBRW= number of reduced velocities per deck 


m = number of decks of input data 


ri 
" 


m 

0.5 + 0.01 IX). (KVBRW). 
De. ), ( F 
A lad 


Machine Components Used 

A number of core storages 

Standard FORTRAN input tape (N1) 
Standard FORTRAN output print tape (N2) 


Standard FORTRAN output punch tape (N3) 
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SECTION VI 
PROGRAM NOTES 


A. Subroutines Used 
BINPU: binary punch routine 
MPUNCH: punch AIC matrix 
MPRINT: prints AIC matrix 
B. Generalized Tapes 


Input, print, and punch tapes are defined as Units 2, 3, and 5, 
respectively; however, these may be altered by placing the desired units 
on symbolic cards HM150007, HM150008, and HM150009. 
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SECTION VIII 
SYMBOLIC LISTING 


Some of the symbols used in the program are defined as follows: 


FORTRAN symbols 


Definition 


A (I,J). [c, .]. also the associated working array 
B (I,J) Lc al also the associated working array 
BR b 

r 
CBAR c 
C (I,J) Korwa also the associated working array 


CH (I,J), CH (I,J + 1) 


DFI, DF2, ..., DF9 


Respective real and imaginary parts of 


the oscillatory AIC Matrix 


Dig: Dy. Das D Digi Dig D 


2’ 3° N-2’ N’ n-]?’ 
Di and Di tp respectively 
DELTAX (1) 4 for segmenti, i=n 
DELVOL (I) AV. for segment i, i=n 
Ix Number of control points 
S s 


LOCAL AREA (I) 


Used in printed output only, see SVAR (below) 


SREF S 

SVAR (I) Sn-1/2 for segment i (i= 1, IX) and 
Sn41/2 for segment IX 

SVMINU (I) Sn-1/2 for segment i 

KVBRW Number of reduced velocities included in 


the data deck 
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FORTRAN symbols Definition 


X (I) Xy-1/2 for segment i (i= 1, IX) and 
41/2 for segment IX 

XKR (J) V/k,, for reduced velocity j, j} = 1, KVBRW 

XN (I) Control point coordinate, segment i 

XMINUS (1) K-1/2 for segment i 

XPLUS (I) Xt1/2 for segment i 


The complete symbolic listing is given on the following pages. 
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